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ABSTRACT 

We perform high resolution ray-tracing simulations to investigate probability distribution functions 
(PDFs) of lensing convergence, shear, and magnification on distant sources up to the redshift of z s = 
20. We pay particular attention to the shot noise effect in TV-body simulations by explicitly showing 
how it affects the variance of the convergence. We show that the convergence and magnification PDFs 
are closely related with each other via the approximate relation /i = (1 — k)~ 2 , which can reproduce the 
behavior of PDFs surprisingly well up to the high magnification tail. The mean convergence measured 
in the source plane is found to be systematically negative, rather than zero as often assumed, and 
is correlated with the convergence variance. We provide simple analytical formulae for the PDFs, 
which reproduce simulated PDFs reasonably well for a wide range of redshifts and smoothing sizes. 
As explicit applications of our ray-tracing simulations, we examine the strong lensing probability and 
the magnification effects on the luminosity functions of distant galaxies and quasars. 
Subject headings: cosmology: theory - gravitational lensing - large-scale structure of universe - meth- 
ods: N-body simulations 
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1. INTRODUCTION 

When light rays from distant sources propagate 
through the inhomogeneous matter distribution in the 
Universe, they are scattered many times by intervening 
clumps of matter. Because of this gravitational lens- 
ing effect (e.g. Schneider, Ehlers & Falco 1992; Schnei- 
der, Kochanek & Wambsganss 2006), a light bundle that 
propagates through the overdense (underdense) region is 
magnified (demagnified) and is also deformed. There- 
fore, observed brightnesses, sizes, and shapes of distant 
objects can in fact differ from those computed assuming 
the homogeneous universe, suggesting that all cosmolog- 
ical observations are subject to gravitational lensing. 

Gravitational lensing effects arc more pronounced for 
more distant sources. Thus it is of growing important in 
modern astronomy where many high redshift objects are 
being discovered. For instance, Bouwens et al. (2010a, b) 
reported candidates of galaxies at z > 7 in deep images 
taken with the Hubble Space Telescope, and derived the 
luminosity function at the redshift (see also Oesch et al. 
2010; Yan et al. 2010). Lehnert et al. (2010) presented 
the spectroscopy of a distant Lyman-a emitting galaxy, 
whose redshift of z = 8.6 turned out to be higher than 
those of previously known most distant objects, including 
a gamma-ray burst at z = 8.2 (Salvaterra et al. 2009; 
Tanvir et al. 2009) and a galaxy at z = 6.96 (lye et al. 
2006). The recent progress in the exploration of the deep 
universe made it possible to derive accurate luminosity 
functions of z-dropout galaxies and Lyman-a emitters at 
z ~ 7 (Ouchi et al. 2009, 2010; Ota et al. 2010). It is 
clear that one has to take proper account of gravitational 
lensing effects for studying intrinsic properties of these 
very distant sources (see, e.g., Wyithe et al. 2011, for a 
recent study of the lensing effects on very high redshift 



galaxies at z > 10). 

Turning the problem around, we can even make use 
of lensing magnification as a 'cosmic telescope' to detect 
and explore faint high-redshift objects that are not acces- 
sible without the help of lensing. Indeed, several highly 
magnified galaxy candidates at up to z ~ 10 have been 
discovered behind massive clusters (Bouwens et al. 2009; 
Bradley et al. 2008; Richard et al. 2006, 2008; Stark et 
al. 2007; Bayliss et al. 2010). Lensing magnification can 
provide a way to locate very high-redshift supernovae 
as well (e.g., Oguri, Suto & Turner 2003; Goobar et al. 
2009; Oguri & Marshall 2010). Upcoming wide-field sur- 
veys such as Hyper Suprime-Cam (HSC; Miyazaki et al. 
2006), the Dark Energy Survey (DES) 1 and Large Syn- 
optic Survey Telescope (LSST; LSST Science Collabora- 
tions et al. 2009) will find more such highly magnified 
distant objects, which serve as good follow-up targets for 
the next generation telescopes such as the Thirty Meter 
Telescope (TMT) 2 and the James Webb Space Telescope 
(JWST) 3 . 

Recently lensing magnification attracts much atten- 
tion in the survey of submillimeter galaxies. Because 
of the steep number counts of submillimeter galaxies, it 
is expected that their observed number counts are sig- 
nificantly modified particularly at the bright end such 
that bright submillimeter sources are almost exclusively 
highly magnified lensing events. This prediction has been 
confirmed by Negrello et al. (2010), who showed that 
bright submillimeter sources from the Herschel Astro- 
physical Terahertz Large Area Survey are indeed domi- 
nated by strong lensing events. Similar enhancements of 
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the number density at the bright end was pointed out by 
van der Burg et al. (2010) for the UV galaxy luminosity 
function at z = 3 — 5 from the Canada-France-Hawaii- 
Telescope Legacy Survey (CFHTLS). 

Furthermore, measurements of distances to any dis- 
tance indicators, which serve as one of the most fun- 
damental ways to constrain cosmological parameters in- 
cluding dark energy, are always affected by gravitational 
lensing. For instance, gravitational lensing not only in- 
duces an external dispersion in distance-redshift relations 
derived from standard candles/sirens such as type la su- 
pernovae, gamma ray bursts, and the gravitational waves 
from neutron star binaries and binary black holes (e.g., 
Hamana & Futamase 2000; Hirata et al. 2010; Shang 
& Haiman 2011), but also may cause systematic biases 
in the derived relations given the non-Gaussian nature 
of the probability distribution of lensing magnification 
(e.g., Sarkar et al. 2008). Kronborg et al. (2010) found 
an evidence of lensing in the distant type-la supernovae 
with a 2cr significance level, from the analysis of the Su- 
pernova Legacy Survey data. The external convergence 
caused by lensing is also known as one of the most sig- 
nificant sources of systematic effects in constraining the 
Hubble constant from time delays between quasar images 
(e.g., Oguri 2007; Suyu et al. 2010) 

In order to address the effects of gravitational lens- 
ing mentioned above, we need an accurate and reliable 
model of the probability distribution function (PDF) of 
gravitational lensing. A powerful way to predict the 
PDF in a given cosmological model is to resort to the 
ray-tracing simulations. In fact ray-tracing simulations 
have been used to study a wide range of gravitational 
lensing effects, including the magnification PDF P(/i) 
(Refsdal 1970; Schneider & Weiss 1988a, b; Jaroszynski 
et al. 1990; Wambsganss et al. 1998; Tomita, Asada 
& Hamana 1999; Tomita, Premadi & Nakamura 1999; 
Barber et al. 2000; Hamana et al. 2000; Wang et al. 
2002; Takada & Hamana 2003; Yoo et al. 2008), the 
convergence PDF P(k) (Jain et al. 2000; Taruya et al. 
2002; Das & Ostriker 2006), the shear PDF P( 7 ) (Bar- 
ber et al. 2000; Jain et al. 2000), the relations among 
the magnification, convergence and shear (Barber et al. 
2000; Hilbert et al. 2011), the strong lensing probabil- 
ity (Wambsganss et al. 1995; Bartelmann et al. 1998; 
Meneghetti et al. 2005; Hilbert et al. 2007), two point 
angular correlation function of the magnification (Takada 
& Hamana 2003), and the cosmic shear (Jain et al. 2000; 
Hamana et al. 2002; Vale & White 2003; White & Vale 
2004; Semboloni et al. 2007; Hilbert et al. 2009; Sato et 
al. 2009, 2011; Sato, Ichiki & Takeuchi 2011). However, 
accuracies of many of previous ray-tracing simulations 
are limited by the resolution. While we have to resolve 
the sub-galactic scales (< 10 kpc) in order to make a 
reliable prediction for the magnification of distant galax- 
ies, only two simulations reached the galactic-scale res- 
olution. One is the simulation by Wambsganss et al., 
who used the small box 5/i _1 Mpc on a side with 256 3 
dark matter particles, the mean particle separation is 
20/i~ 1 kpc and the softening length is 10/i _1 kpc (Cen et 
al. 1994; Wambsganss et al. 1995, 1997, 1998). They cal- 
culated the strong lensing probability, the distribution of 
the magnification, and the lensing effects on the determi- 
nation of the deceleration parameter qo using type-la su- 
pernovae. The other one is the work of Hilbert et al., who 



used the Millennium Simulation (Springel et al. 2005) in 
which the box is 500/i _1 Mpc on a side with 2160 3 parti- 
cles, the mean particle separation is 230/i~ 1 kpc and the 
softening length is 5k" 1 kpc (Hilbert et al. 2007, 2008, 
2009). They discussed many of the issues such as P(n), 
P(k), the strong lensing probability and the effects of 
baryon. 

An alternative method to predict lensing PDFs is to 
use the so-called halo model (e.g., Jain &; Lima 2011; 
Lima et al. 2010a,b; Kainulainen & Marra 2009,2011, 
for a recent example). While the halo model can make 
accurate predictions of PDFs at the high magnification 
tail where contributions from single halos dominate, pre- 
dicting the PDFs near the peaks with the halo model 
is quite challenging. Thus, cross-checking of the halo 
model with high-resolution ray-tracing simulations are 
crucial for validating and possibly improving the accu- 
racy of the halo model predictions. 

In this paper, we perform high resolution ray-tracing 
simulations to study lensing PDFs in great detail. The 
box size of 50h~ 1 Mpc with 1024 3 particles, the mean 
particle separation of 50/i _1 kpc, and the softening length 
of 2/i~ 1 kpc indicates that our simulations represent, 
to our knowledge, the highest resolution ray-tracing 
simulation conducted to date for studying cosmological 
lensing effects. We run three sets of simulations with 
256 3 , 512 3 and 1024 3 particles to check the numerical 
convergence of our simulation results, and also compare 
the two box sizes of 50/i _1 Mpc and 100/i _1 Mpc to 
check the effects of density fluctuation larger than the 
box size. We consider the gravitational evolution of 
the dark matter particles only, and do not include the 
effect of baryon cooling which would enhance the strong 
lensing probability (e.g., Hilbert et al. 2008). In this 
paper, we present comprehensive analysis of PDFs of 
convergence, shear, and magnification, with a particular 
emphasis on the relation between these three lensing 
quantities. Our exploration of the PDF up to very high 
redshift of z s ~ 20 enables the immediate application 
of our results for studying high-redshift (z > 7) sources 
that are recently discovered and also predicting even 
more distant sources that will be discovered in the 
future. Our numerical results of two-dimensional maps 
of lensing fields and PDFs are made publicly available at 
http://cosmo.phys.hirosaki-u.ac.jp/takahasi/raytracing/. 

The structure of the present paper is as follows. In 
Section [21 we describe the details of iV-body simulations 
and our ray-tracing simulations and method to obtain 
two-dimensional maps of the convergence, the shear and 
the magnification. In Section [3l we show the PDFs of 
the convergence, the shear and the magnification and 
provide simple analytical formulae of these PDFs which 
well reproduce PDFs obtained from ray-tracing simula- 
tions over a wide range of redshifts and smoothing scales. 
In Section |4j we examine the strong lensing probabil- 
ity by categorizing lensing mapping into tree types. In 
Section [5j we discuss magnification effects on luminos- 
ity functions of distant sources. Section [6] is devoted to 
summary and discussion. 

Throughout the present paper, we adopt the standard 
ACDM model with matter density Vt m = 0.274, baryon 
density = 0.046, dark energy density fi w = 0.726, 
spectral index n s = 0.96, amplitude of fluctuations 
erg = 0.812, and expansion rate at the present time 
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Fig. 1. — A schematic picture of our ray-tracing simulation. The 
horizontal axis is the comoving distance r from the observer. The 
vertical thick lines denote the positions of the lens planes, which 
are located at r = L X (i-f 1/2) with i = 0, 1, 2, ■■. The source planes 
are located at r = L X i. Light rays are emitted from the observer 
and are scatted at the lens planes before reaching the source plane. 
The field of view is 1 X ldeg 2 . 

Ho = 70.5km s" 1 Mpc" 1 , consistent with the WMAP 
5-year results (Komatsu et al. 2009). 

2. NUMERICAL SIMULATIONS 

Figure [T] shows a schematic picture of our ray-tracing 
simulation. The horizontal axis is the comoving dis- 
tance r from the observer, and thick vertical lines are 
the lens planes. The lens planes are placed at an equal 
distance intervals of L, r = L x (i + 1/2) with an inte- 
ger i = 0,1,2,- • •. We place the lens planes from the 
observer to the highest redshift of z s = 20 with the in- 
terval L = 50/i _1 Mpc and 100/i _1 Mpc. Light rays are 
emitted from the observer and are deflected at each lens 
plane before reaching the source plane. Source planes are 
placed in between the lens planes, i.e., at r — L x i. In 
our simulation, the field of view is 1 x ldeg 2 . We impose 
the periodic boundary condition on the lens planes. 

2.1. N-body Simulations 

We run the TV-body simulations on the cubic box, and 
then project the particle positions to the two dimen- 
sions, in order to obtain the particle distribution and the 
gravitational potential on the lens planes. We use the 
numerical simulation code Gadget2 (Springel, Yoshida 
& White 2001; Springel 2005). We generate the ini- 
tial conditions based on the second-order Lagrangian 
perturbation theory (2LPT; Crocce, Pueblas & Scocci- 
marro 2006; Nishimichi et al. 2009) with the initial lin- 
ear power spectrum calculated by Code for Anisotropies 
in the Microwave Background (CAMB; Lewis, Challi- 
nor & Lasenby 2000) 4 . We use a parallelised 2LPT 
code which is kindly provided by Takahiro Nishimichi 
(Valageas & Nishimichi 2011) to run large cosmologi- 
cal iV-body simulations with initial conditions based on 
2LPT. We dump the outputs (the particle positions) at 
the redshifts corresponding to the positions of the lens 
planes r = L x (i + 1/2), shown in Figure [T] The size 
of the simulation box is L on a side. We have checked 
that the matter power spectra of our iV-body simulations 
agree with the results of the higher resolution simulation, 
in which we used the finer simulation parameters of the 
time step, the force calculation, etc., within 2(10)% for 
k < 20(80)/i/Mpc. 

4 see also http://camb.info/ 



Table Q] lists four models of our simulations. The ta- 
ble shows the simulation box on a side L, the number of 
particles iV 3 , the particle mass m p , the softening length 
r s , the initial redshift z- m , and the number of lens planes 
N\ cns up to z s — 20. The softening length is fixed to be 
4% of the mean particle separations for all the models. 
The first three models are L = 50/i _1 Mpc with differ- 
ent number of particles, 256 3 , 512 3 , and 1024 3 , which we 
name S256, S512, and S1024, respectively. We use these 
three models to check the numerical convergence of our 
simulation results at small scales, particularly because 
the lensing magnification is known to be very sensitive 
to the mass resolution. As we will show later, the poor 
mass resolution simulations are indeed affected signifi- 
cant by the shot noise. The model L512, L = 100/i _1 Mpc 
with Np = 512 3 , is the same resolution as S256 but the 
simulation volume is eight times larger than S256. We 
use this model to check the effect of density fluctuation 
larger than the box size (see Appendix A). For S256, 
S512 and L512, we prepare the different realizations at 
each lens plane to reduce the sample variance 5 . How- 
ever, for S1024, we perform only four independent real- 
izations because of the limited computer resources, and 
we repeatedly use the outputs at different redshifts to 
construct a light-cone output. 

2.2. Ray-tracing Simulations 

We briefly explain the procedure to trace light rays 
through TV-body data and obtain the maps of the lensing 
fields on the source plane. We use the code RAYTRIX 
(Hamana & Mellier 2001) which follows the standard 
multiple lens plane algorithm. In the standard multi- 
ple lens plane algorithm, the distance between observer 
and source galaxies is divided into several intervals. In 
our case, as shown in Figure [T] we adopt a fixed in- 
terval whose value is the same as simulation box L on 
a side. Particle positions are projected onto two di- 
mensional lens planes (xy, yz ,zx planes) every L. Us- 
ing Triangular-Shaped Cloud method (Hockney & East- 
wood 1988), we assign the particles onto Ng grids in lens 
planes, then compute the projected density contrast at 
each plane. We test the convergence of our simulation by 
varying resolution from TV 2 = 256 2 to 16384 2 . The two- 
dimensional gravitational potential is solved via Poisson 
equation using Fast Fourier Transform. Finally, two di- 
mensional sky maps of the convergence, the shear, the 
magnification, and the angular positions of light rays are 
obtained by solving the evolution equation of Jacobian 
matrix along the light-ray path which is obtained by solv- 
ing the multiple lens equation. 

We prepare 20 realizations by randomly choosing the 
projecting direction and shifting the two dimensional po- 
sitions. In each realization, we emit 2048 2 light-rays, 
leading to 8 x 10 7 rays in total. 

We note that our work is complementary to the pre- 
vious work by Sato et al. (2009, 2011) who conducted 
ray-tracing simulations using the similar technique as de- 
scribed above. Their interests lie in accurate predictions 
for the covariance matrices of the cosmic shear power 
spectrum and correlation function, which are required 

5 We run the TV-body simulation for the i-th lens plane at r = 
L X (i — 1/2) only from the initial redshift to the redshift of the 
lens plane. 
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TABLE 1 



Models 


L (ft-J-Mpc) 


TV 3 
p 


m p (h- 1 M Q ) 


r s (h i kpc) 


Zin 


-^lens 


S256 


50 


256 3 


5.7 x 10 s 


8 


80 


158 


S512 


50 


512 3 


7.1 x 10 7 


4 


90 


158 


S1024 


50 


1024 3 


8.9 x 10 6 


2 


100 


158 


L512 


100 


512 3 


5.7 x 10 s 


8 


80 


79 



Note. — Models of our ray-shooting simulations. The cubic simulation 
box on a side L, the number of particles JV 3 , the particle mass m p , the 
softening length r s , the initial redshift Zi n , and the number of lens planes 

-^lens U P to Z s = 20. 



for extracting cosmological information from future wide 
field optical imaging surveys (e.g., HSC survey, DES and 
LSST), whereas our purpose is to predict the lensing 
PDFs for distant sources, for which high resolution sim- 
ulations are crucial. 

3. RESULTS 
3.1. Variance of the Convergence 

The surface area of the ray bundle in the source plane 
is smaller (larger) than that in the image plane due to the 
cosmic magnification (demagnification). As a result, the 
probability distributions of the convergence, shear and 
magnification evaluated in the image plane differs from 
those in the source plane. Throughout this paper, we 
will show the PDFs defined in the source plane, because 
our interest lies in the predictions of lensing effects on 
the distant sources, for which source-plane lensing PDFs 
are more relevant. Since in our ray-tracing simulation 
we emit the light rays homogeneously from the observer 
to the first lens plane, we can easily derive the source 
plane PDFs by adding a weight of the magnification in 
computing the PDF from simulations (see below). 

In this section, we show our numerical results of the 
variance of the convergence k in the source plane : 

Here Kj and /Xj denote the convergence and magnification 
for j-th light ray, respectively. The summation runs over 
all the light rays of 8 x 10 7 . The factor 1/fij originates 
from the ratio of the area in the image and source planes. 

Figure [2] shows the root-mean-square of the conver- 
gence, (k 2 ) 1 / 2 . The horizontal axis is the grid size of 
the two dimensional gravitational potential in the lens 
planes. Here the grid size is r gr id = L/N g with the num- 
ber of grid Ng of the gravitational potential, and our 
finest resolution is r gr jd = 3ft. _1 kpc (corresponding to 
Ng = 16384 2 ). The symbols are our simulation results, 
the cross x, the circle • and the plus +, correspond to 
the various number of particles, iVj = 256 3 (model S256), 
512 3 (S512) and 1024 3 (S1024), respectively. The source 
rcdshifts are z s — 1, 2, 3, 5, 10 and 20. 

In the weak lensing approximation, the variance of the 
convergence is given by (e.g., Bartelmann & Schneider 
2001), 



< K 2 } = ±Hfa 2 m 
an 



dz 



(i + zy 



r(z)r(z,z s ) 



r{z s ) 
dkkP(k,z)(2) 



where r(z, z s ) is the comoving distance from z to z s and 
P(k, z) is the matter power spectrum as a function of 
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Fig. 2. — Root-mean-square of the convergence, (k 2 ) 1 / 2 , as 
a function of the smoothing scales r gr i<j at the source redshifts 
of z s = 1,2,3,5,10 and 20. Different symbols denote simulation 



results with different mass resolutions, the cross X for N„ 



256 3 



the circle • for 512 3 and the plus + for 1024 3 , respectively. The 
dashed lines are the weak lensing prediction given in Equation J2} . 



the wavenumber k and z. We calculate the matter power 
spectrum P(k, z) in Equation (J2J) directly using our TV- 
body simulation data at the redshifts of the lens planes 6 . 
The predictions are shown by horizontal dashed lines in 
Figure O We do not use the theoretical fitting formula 
of the halo-fit model (Smith et al. 2003) to calculate 
the non-linear power spectrum, because it has recently 
been reported that the halo fit underestimates the power 
spectrum at small scales, k > 0.1/i/Mpc (e.g. White & 
Vale 2004; Sato et al. 2009; Heitmann et al. 2010). For 
reference, the variance in Equation ([2} becomes smaller 
by a few ten percents if we use the halo-fit model to 
compute the non-linear matter power spectrum P(fc, z) 
in Equation (f2]). 

As shown in Figure [2 the simulation results decrease 
for larger r gr id- This is because the density fluctuations 
smaller than r gr jd are smeared out when assigning the 
particles onto the grids. On the other hand, for smaller 
r gridj our simulation results are larger than the theo- 
retical prediction (the dashed lines) especially for the 
higher redshifts and for the lower mass resolution sim- 
ulations, which we ascribe to the shot noise. For higher 

6 To calculate the density fluctuations in the cubic box, we as- 
sign the particles onto a 1280 3 grid using the cloud-in-cell method. 
Then we perform the Fourier transform to calculate the power spec- 
trum. 
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Fig. 3. — Same as previous figure l2l for the model S256 (i.e., the 
lowest resolution run), but including the effects of the smoothing 
and the shot noise in the theoretical model. The solid curves are 
the theoretical model including effects of both the shot noise and 
the smoothing in Equation (3). The dotted curves are predictions 
with the smoothing effect only. 



source redshift, the smaller density fluctuations gener- 
ate the convergence k in Equation (the peak position 
of integrand in Equation @, k 2 P(k,z), shifts larger k 
for higher z). Hence, in order to study the lensing of 
high redshift sources using the ray-tracing simulations, 
we need sufficient number of particles to reduce the shot 
noise. 

To understand the effects of the shot noise and the 
smoothing, we rewrite Equation ([2]) by taking these ef- 
fects into account: 



<- 2 > 



dz (~\ , \ 2 



r(z)r(z, z s 



dkk 



P(k,z) 



r(z s ) 
e -(fcAc Ut ) 2 j(3) 



where 1/n is the shot noise term (n is the number den- 
sity of the iV-body particles defined as n = N 3 /L 3 ) and 
e -(fe/feout) j s t ne Gaussian smoothing term with a cutoff 
wave number /c cut — ax 2tt /r ST id i a is a constant order 
of unity; we adopt a = 0.2). In Figure [3J the theoretical 
model given by Equation © is compared with the lowest 
resolution ray-tracing results. We find that this model 
with the shot noise and smoothing is in good agreement 
with the simulation results for both the small and the 
large r gr id- For comparison, dotted curves show the pre- 
dictions when only the smoothing term in Equation ([3]) 
is included. The differences between the solid and dot- 
ted curves clearly show the effect of the shot noise at 
small r gr id. Furthermore, our calculation provides a way 
to infer the critical grid scale from which the shot noise 
dominates. We note that we plot the low resolution sim- 
ulation results with N 3 = 256 3 to show the shot noise 
effect clearly. As shown in Figure [2J in the highest res- 
olution model S1024, the shot noise is negligible even at 
redshift z s ~ 20 and the smallest grid size of 3/i -1 kpc. 
In what follows, we use this highest resolution model to 




1+K 

Fig. 4. — Convergence PDF as a function of 1 + re for distant 
sources at z s = 1 — 20. We show the results of our highest resolu- 
tion simulation (model S1024, the smoothing scale 3fe _1 kpc). The 
arrows indicate the minimum convergence for the empty beam. 



compute the PDFs for different redshifts and grid scales. 

3.2. Convergence, Magnification and Shear PDF 

In this section, we show our simulation results of the 
probability distribution of the convergence, magnifica- 
tion and shear up to z s = 20. We will show the PDFs for 
the highest resolution simulation of model S1024 with 
the smoothing length 3ft. _1 kpc. We have checked that 
the PDFs are consistent with those computed from lower 
resolution simulations, as long as the grid size is large 
enough for the shot noise not to dominate. 

As noted in Section 13.11 the PDF in the source plane 
dPs I dx is different from that in the image plane dPj j dx 
by a factor of the magnification: 



dPsix) 
dx 



1 dPr(x) 
M 



dx 



(4) 



where x denotes the convergence n, shear 7, and magni- 
fication //. 

3.2.1. Convergence PDF 

Figure 2] shows the convergence PDF as a function of 
1 + not k itself, at redshifts z s = 1 — 20. For higher 
redshifts, the distribution becomes broader and its peak 
moves to lower value. These features are consistent with 
previous works (e.g., Jain et al. 2000). With increas- 
ing the smoothing scale r gr id, the distribution becomes 
narrower as expected from the variance in Figure [2j As 
shown in the Figure, there are small dumps at 1 + k ~ 2 
because of the multiple images form for n > 1. Given 
the way we compute the source-plane PDFs, we count 
the images more than once for a single source. We will 
discuss how to treat the multiple images in Section |4j 

The convergence has a minimum value when the light 
ray propagates through the empty region (so called the 
empty beam, e.g. Jain et al. 2000). The convergence for 
the empty beam is given by, 



''empty 



--i?oOr. 



im { ] r( Zs ) • (5) 
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Fig. 5. — The negative mean of the convergence as a function 
of the variance of the convergence. The dots are our simulation 
results for various smoothing scales rg r i<j = 3 - 200/1" x kpc and 
source redshifts z s = 1 — 20. The solid line is a simple theoretical 
prediction, (ft) = — 2(ft 2 ) (see text for details). The mean conver- 
gence has a strong correlation with its variance, which is explained 
well by the simple relation above. 



The arrows in Figure [4] show the minimum convergence 
of the empty beam, 1 + K cmp t y . As shown in the Figure, 
for a low source redshift (e.g. z s = 1), the minimum 
convergence of the empty beam (the arrow) is consistent 
with the minimum value of the PDF. However, for higher 
redshifts, the convergence of the empty beam becomes 
significantly lower than the minimum values of the PDFs. 
This is because the matter distribution becomes more 
homogeneous for higher redshift universe and hence there 
is little chance to propagate through the empty region. 

In order to check our numerical simulations, we eval- 
uate the mean of the magnification (/i) and the con- 
vergence (k). We find that the mean magnification is 
(fi) = 1 within small scatters less than 1% for z s = 1 — 20. 
Hence, by averaging many light rays, we recover the 
filled-beam distance assuming the homogeneous mass 
distribution in the standard cosmological model. On the 
other hand, we find that the mean convergences (k) in 
ray-tracing simulations are not zero but systematically 
smaller than zero. Figure [5] shows our numerical results 
of the mean convergence (k) as a function of the variance 
(k 2 ). We show the results for various smoothing scales 
r gr jd = 3 — 200/i~ 1 kpc and source redshifts z s = 1 — 20. 
As shown in the Figure, (k) is systematically negative, 
which is contradictory to the naive expectation in the 
weak lensing limit, (k) = 0. We find that the mean is 
well correlated with the variance, which can be under- 
stood as follow. Assuming k is small, we can expand the 
magnification as, 



1 



(1-k) 2 - 7 2 ' 
~1 + 2k + 3k 2 + 7 2 . (6) 

Then, by taking the mean of the above equation, we 



obtain 

1 = <M) - 
from which we have (k) 



1 



-2(«) - 



(7) 



K"). The solid line in Fig- 
ure 03 corresponds to the theoretical prediction in Equa- 
tion ||7J). We find that the simulation results are ex- 
plained well by this simple relation. The small devia- 
tion from the solid line can be due to the higher order 
moments of the convergence and the shear neglected in 
Equation {BJ. Our simulation indicates that the light 
rays pass through underdense region on an average. This 
area magnification effect on the average value of the con- 
vergence was also derived by a perturbative approach 
(Hamana 2001). A possible explanation for this is that, 
when a light ray passes near a massive object, its paths 
is defected such that the closest distance to the massive 
object becomes larger and hence the ray passes through 
the less dense region. The negative mean of the conver- 
gence, for instance, can have an impact on the statistics 
of time delays (Oguri 2007). 

We note that the above results depend on which plane 
(image or source plane) we use to calculate the mean, be- 
cause the higher magnification events are more weighted 
in the image plane. We find that the mean conver- 
gence in the image plane is zero within small scatters 
less than 10~ 3 . This is because the mean in the image 
plane is different by a factor of fi from Equation Q, 
and hence the mean convergence in the image plane is 
{ft)imagc = — ((1 + 2k)k) = from Equation 0- 

Similarly, the mean magnification in the image plane is 
found to be systematically positive, (/x)i mage = (fi 2 ) > 0. 
In addition, our results are derived under an ideal sit- 
uation in which light sources are distributed uniformly 
in a source plane and all the sources are supposed to be 
observed. In general, only sources with luminosity above 
a certain threshold are observed. In that case, the mean 
magnification can be non-zero with either sign because 
of the so-called magnification bias (e.g. Schneider et al 
1992 and see section 5.2). Hence, we have to use the 
appropriate statistics depending on situations. 

Finally, we compare our numerical results of the con- 
vergence PDF with previous works. Taruya et al. (2002) 
conducted the ray-tracing simulations to investigate the 
statistical properties of weak-lensing field. They showed 
that the convergence PDF is well described by the log- 
normal distribution, which reflects the fact that the one- 
point distribution function of matter density field is well 
described by the lognormal model (Kayo et al. 2001). 
Later, Das & Ostriker (2006) calculated the probability 
distribution of the projected surface mass density using 
cosmological iV-body simulations to show that the PDF 
is well fitted by the modified lognormal distribution: 
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= N K exp 
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1 + K/ |K omp ty| J K + |K cm pty| 

where N K is the normalization. By setting A K = and 
u K = (k 2 ), Equation © reduces to the lognormal model 
of Taruya et al. (2002). The two parameters A K and uj k 
are determined by using the following conditions, 
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Fig. 6. — The convergence PDF for various smoothing scales 
'"grid an d redshifts z s . Black dots are our simulation results, while 
red solid curves are the modified log-normal fitting model in Das 
& Ostriker (2006). 
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We use our simulation results of (n 2 ) in the above equa- 
tions. Note that the first condition of Equation © dif- 
fers from the original condition used in Das & Ostriker 
(2006), who adopted (k) = 0. We modify this condition 
to the form above based on our finding shown in Figure [5] 
and Equation 0. 

Figure [5] is the comparison of the convergence PDF 
between our simulation results and the model of Das & 
Ostriker (2006). Black dots are our simulation results, 
while red curves are their fitting function. The nine pan- 
els are for various smoothing scales of r gr id = 3hr 1 kpc 
(left), 12/i~ 1 kpc (middle) and 48/i _1 kpc (right) and for 
various source redshifts of z s = 1 (top), 3 (middle) and 10 
(bottom). As seen in the Figure, the fitting model agrees 
well with our results in all the panels, especially near the 
peak of the distribution. However, for high convergence 
tail, the model slightly underestimates the PDFs 

3.3. Magnification PDF 

Figure [7] shows the magnification PDF of the distant 
sources at z s = 1 — 20. We use the highest resolu- 
tion model S1024 (see Table[T| with the smoothing scale 
3/i _1 kpc. As clearly seen in the Figure, for more distant 
sources, the peak of PDF moves to the fainter magnifi- 
cation (less than 1) and its distribution becomes broader 
as in Figure These features are consistent with the 
previous works (e.g., Wambsganss et al. 1998; Wang et 
al. 2002; Hilbert et al. 2007). For example, at the high- 
est redshift z s = 20, 10%(1%) sources are magnified by 
over 30%(100%). For high magnification (fx > 10), the 
simulation results are asymptotically proportional to ji~ 2 
which is consistent with the analytical expectation (e.g., 
Schneider et al. 1992). We note that there are small in- 
creases of PDFs at very low magnification (/i < 1). These 
correspond to fainter images of strongly lensed multiple 




Fig. 7. — Magnification PDFs of the distant sources at z s = 
1 — 20. We show the results of our highest resolution simulation 
(model S1024, the smoothing scale 3/i _1 kpc). 
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Fig. 8. — The two dimensional distribution of the convergence 
and the magnification. The solid curve is jtt = (1 — t)~ 2 and the 
dashed line is fi = 1 + 2re (weak lens approximation). 



images. Again, we discuss strong lensing events in more 
detail in Section |4j 

In previous works, the convergence and magnification 
PDFs have been presented separately. Here we inves- 
tigate the relation of these two PDFs. First we check 
the correlation of the convergence k and magnification 
H, which is shown in Figure [8] The vertical axis is fi, 
while the horizontal axis is k, for redshifts z s — 1 — 20. 
The Figure indicates that convergence and magnifica- 
tion are highly correlated with each other, which is con- 
sistent with Barber et al. (2000) and Hilbert et al. 
(2011). We find that the correlation is well explained by 
/it = (1 — k)~ 2 (solid curve), which comes from the defini- 
tion of the magnification with the shear term neglected. 
On the other hand, the relation \i = 1 + 2k (dashed line) , 
which holds in the weak lensing approximation, cannot 



Takahashi et al. 



tion PDF as 



=3kpc/h 12kpc/h 



48kpc/h 



CL 

T3 




Fig. 9. — The magnification PDFs for various smoothing scales 
r grid and redshifts z s . Black dots are our simulation results. Or- 
ange dots denote the magnification PDFs converted from the con- 
vergence PDFs in our simulations via the relation fi = (1 — k) -2 . 
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Fig. 10. — The magnification PDF for various smoothing scales 
r grid an d redshifts z s . Black dots shows our simulation results, 
green dotted curves are the fitting model in Wang et al. (2002), 
red solid curves are the fitting function in Das & Ostriker (2006) 
converted via the relation fi = (1 — ft) - 2 , orange dashed curves are 
same as the red solid curves but including the high magnification 
tail. 



explain the correlation in our simulations very well. 

Given the tight correlation, we translate the conver- 
gence PDF to the magnification PDF using the relation 
[i = (1 — k)~ 2 . Specifically, we compute the magnifica- 



dP s 



(1 



dP x 



dn 



(10) 



The result is shown in Figure The black dots are the 
original magnification PDFs, while the orange dots are 
those converted from the convergence PDFs in our sim- 
ulations using the relation /i = (1 — k)~ 2 . Note that we 
use the simulated convergence PDFs only for k < 1, as 
the conversion equation clearly breaks down at k = 1. 
We find that these two PDFs agree surprisingly well, up 
to very high magnification tails of /x ~ 10. Hence, once 
either the convergence or magnification PDF is available, 
we can easily obtain the other PDF by the transforma- 
tion given above. 

We now compare our numerical results of magnification 
PDF with previous works. Wang et al. (2002) proposed 
the fitting function of the magnification PDF using sim- 
ulation results of Wambsganss et al. (1997), Barber et 
al. (2000) and Munshi & Jain (2000). Their formula 
depends only on the variance of the convergence, and in- 
dependent of cosmological models and redshifts. Their 
fitting function is a stretched Gaussian distribution. In 
Figure 1 101 green dotted curves show their model, while 
black dots are our simulation results. We find that their 
model agree well only near the peak of the distribution. 
In fact, as they noticed, their model can be used only 
for near the peak fi ~ 1. They used the approximation 
of fj, = 1 + 2k and their model depends on the variable 
Mmin = 1 + 2K emp t y . However, for high redshift z s > 4, 
this quantity ^ m i n becomes negative, leading to the break 
down of the formula (hence we do not plot the green dot- 
ted curves for z s = 10 in Figure fTO]) . 

Since the modified lognormal model in Das & Ostriker 
(2006) well reproduces the convergence PDF, we convert 
this model to the magnification PDF using Equation fj 10|) 
to see if the model reproduces the simulated magnifi- 
cation PDFs. We find that this model (red curves in 
Figure 1 10[) agrees reasonably well with our simulation 
results, although the model cannot reproduce the large 
magnification tail, presumably reflecting the slight un- 
derestimate of the convergence PDFs at the high conver- 
gence tail. Since we know that the magnification PDF 
behaves as dPg/d/i cx [i~ 3 at the tail, we propose a sim- 
ple phenomenological model which takes account of the 
tail behavior, 



dP s 
d/j, 



(1 - Ky dP s 

2 dn 



+ Q(fi — 1) exp 

, (1-k ) 3 dP s 
< 2 dn 
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40- ir 
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where O is the step function and /iq = (1 — kq)~ 2 . We 
set no = 3 as a reasonable choice to fit our simulation 
results for a wide range of parameters. In right hand 
side of Equation (fTTj) . the convergence PDF dP s /dn is 
evaluated by using the modified lognormal model of Das 
& Ostriker (2006). This model, shown by the orange 
dashed curves in Figure 1101 reproduces the simulation 
results well both near peaks and at high magnification 
tails. We note that a more robust and accurate model 
of the magnification PDF can be obtained by computing 
the tail distribution separately using the halo model to 
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Fig. 11. — The shear PDF of distant sources at z s = 1-20. The 
results are obtained from our highest resolution simulation (model 
S1024, the smoothing scale 3h~ 1 kpc). 
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Fig. 12.— The shear PDF for various smoothing scales r grid 
and redshifts z s . Black dots are our sim ulati on results, while red 
curves are our fitting formula in Equation (1 1 2t . Our fitting formula 
agrees with the data very well. 



combine it with the peak distribution computed from the 
modified long-normal convergence PDF. 

3.4. Shear PDF 

Figure [TT] shows the shear PDF as a function of the 
shear amplitude 7 = (7^ + 7! ) 1//2 . With increasing the 
redshift, the peak of the PDF becomes larger while its 
shape is not changed very much. The shear PDF is pro- 
portional to 7 2 for the small shear amplitude (7 -C 1). 
This behavior has been explained by an analytical calcu- 
lation for light rays propagating through randomly dis- 
tributed point-mass lenses (e.g., Schneider et al. 1992, 
Sec. 11.2). The tidal field of the lenses generates the 



shear, while the convergence is exactly zero because the 
light propagates the empty space between the point-mass 
lens particles. The analytical PDF is dP s /d\n^ oc 7 2 for 
7 <C k* and oc 7 _1 for where is the conver- 

gence for the surface density of the uniformly distributed 
point-mass lenses. Our simulation result is consistent 
with the simple theoretical model only for small 7. 

Here we present a fitting formula for the shear PDF. 
We empirically derives the following fitting formula for 
the PDF of 7, 



G?7 



N. 
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exp 



{ln(l + A 77 )} 2 
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(12) 



where -ZV 7 is a normalization. In order to determine the 
parameters A 1 and lj 7 , we use 18 shear PDFs for three 
different smoothing scales r gr jd = 3, 12, 48/i _1 Mpc with 
six redshifts z s = 1, 2, 3, 5, 10 and 20. We find the best- 
fitting values of A 7 and w 7 are well correlated with the 
shear variance a 7 (= (7 2 ) 1 / 2 ) as 

0.010 



(0.00075 + ct 2 ) 
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tcu = 1.66 a. 



1.32 



(13) 



Figure [T^] shows the above fitting formula (red curves) 
and the simulation results (black dots) for various 
smoothing scales r gr id and redshifts z s . As shown in 
the Figure, our fitting formula agrees with the simula- 
tion results quite well for a wide range of redshifts and 
the smoothing scales. 

Finally, we investigate the correlations between shear 
and convergence 

the quantities x and y is defined as, 



A correlation coefficient r xy between 



((x-(x))(y-(y))) 

((x-( a; )) 2 ) 1 /2((y„(y))2 ) l/2' 



(14) 



where x,y represents k, 71,2 and 7. 

We find a positive correlation between k and 7, while 
we find no correlations among 7 1; 72, and k. Figure [TBI 
shows the correlation between k and 7. The vertical axis 
is the correlation coefficient r Ki7 defined in Equation (|14[) . 
and the horizontal axis is the source redshift. The Figure 
shows the strong correlation of 0.1 — 0.5 with a clear ten- 
dency that the correlation is stronger for lower redshift 
sources. 

The positive correlation between k and 7 is expected 
at least when a single lensing event near a halo is domi- 
nant. On the other hand, for multiple lensing by differ- 
ent halos, the total convergence is simply given by the 
sum of each event, whereas the shear can be cancelled 
out depending on the relative direction of the shear in 
each event. This implies that the correlation becomes 
weaker when there are more intervening halos that can 
potentially contribute to lensing. This simple considera- 
tion appears to be consistent with the simulated results 
above, as higher redshift sources have more intervening 
halos that can significantly affect the light propagation. 

4. STRONG LENSING PROBABILITY 

In this section, we calculate the strong lensing proba- 
bility and investigate the effects of the multiple images 
on the PDFs. 

The light-ray path from the source to the observer is 
the stationary point of the time delay r (or the Fermat 
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Fig. 13. — Correlation between 7 and Ki clS EL function of the 
source redshift. The vertical axis is the correlation coefficient de- 
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Fig. 14. — Strong lensing probability for type II (left panel) and 
type III (right panel). 



potential), i.e., dr/dOi = 0, due to the Fermat's prin- 
ciple. The lensing mapping from the image position to 
the source position is characterized using the Jacobian 
matrix which is the Hessian of r, Aij = d 2 r / d9id9j. Us- 
ing the convergence n and the shear 71,2, the Jacobian 
matrix is expressed as (e.g., Schneider et al. 1992) 



A = 



1 - K - 71 -72 
-72 1 - K + 7x 

The mapping is categorized into three types: the min- 
imum, the maximum, and the saddle point of r. 



(15) 



The 



minimum point is called type I, the saddle point is type 
II, and the maximum point is type III (e.g., Schneider et 
al. 1992): 

type I : detA > k tiA > 0, 

i.e., |1 — k| > 7 & k < 1 



type II : detA < 0, i.e., |1 - k\ < 7 (16) 
type III : detA > & txA < 0, 

i.e., |1 — k|>7&k>1 

Type II and III correspond to multiple images produced 
by strong lensing. The strong lensing probability is de- 
fined as the ratio of the number of light rays in theses 
two types, 

p Ei( 1 M)typcII,III f 7 , 

2-(iWM» Mtypes 

The factor 1/fii is due to the probability in the source 
plane. Figure [TH shows the strong lensing probability of 
type II (left panel) and type III (right panel) as a function 
of r gr id. As shown in the Figure, as the smoothing scale 
decreases, the probability increases. This is because the 
smoothing effect smears out the central cusp of the dark 
halo. The probability for type II is higher than that for 
type III. 

Figure [15] shows the contribution of each type to the 
PDFs at z s — 5. The left panel is the magnification 
PDF, the middle panel is the convergence PDF, and the 
right panel is the shear PDF. The green circles are type 
I, the blue crosses are type II, the red pluses are type 
III, and the black dots are the sum of them. As clearly 
seen in the Figure, almost all light rays are type I. With 
increasing the source redshift, the fractions of type II 
and III increase. In the magnification PDF, the small 
increase seen at /z < 0.6 is due to type III. In the high 
magnification limit, type I and II have same probability. 
In the convergence PDF, the small knot at 1 + n ~ 2 is 
type II, and the second knot at 1 + k ~ 3 is type III. In 
the shear PDF, type II and III appear in the high shear 
limit. The resolution of our simulation is high enough to 
see these strong lensing features clearly in the PDFs. 

In order to understand these three types more clearly, 
we consider the strong lensing by a NFW halo (Navarro, 
Frenk & White 1997). For the NFW profile, the caustic 
is the circle on the source plane. If the source position is 
outside the caustic, only one image of type I is formed. 
When the source position crosses the caustic inward, two 
new images are formed. One is formed the opposite side 
of the lens center, which is type II, and the other one is 
a faint image formed close to the center, which is type 
III. Hence, as shown in Figure [TSJ the convergence of the 
type III is very high k > 1 while the magnification is 
small. When the source is close to the lens center, the 
two bright images of type I and II are strongly magnified 
with the similar magnification, which can explain the 
high magnification tail in the left panel of Figure [15] 

5. MAGNIFICATION EFFECTS ON LUMINOSITY 
FUNCTIONS OF DISTANT SOURCES 

5.1. Effects of Source Size on Magnification PDF 

So far we have discussed lensing effects for point 
sources. In this section, we investigate the effects of 
the finite source size on the magnification PDF. The fi- 
nite size smooths the magnification profile on the source 
plane, and smears out some high magnification events 
(/i> 1). The smooth magnification \x of the extended 
source with surface brightness profile l(r\) is generally 
given by, 
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Fig. 15. — Magnification PDF (left panel), convergence PDF (middle panel) and shear PDF (right panel) for z s = 5 for different types 
of lensing mapping: green circles are type I, blue crosses are type II, red pluses are type III, and black dots are the sum of them. 




Fig. 16. — The magnification PDF for the extended source 
at z B = 5. The red (blue) symbols are the results for the source 
radius of 10 (50) h~ J kpc, while the black symbols are for the point 
sources. Here, we do not include the multiple images. 



where r\ — (771,772) is the two-dimensional vector in the 
source plane. For simplicity, we assume that the surface 
brightness is a circular top-hat model with the radius Rs, 
I{r\) = l/(vri?|) for r\ < R s and I(rj) = for 77 > R s . 
Note that the distances rj,Rs are the comoving scale. 
We adopt Rs = 10/i _1 kpc as a typical size of galaxies. 

In order to resolve the sub- galactic scales (< 
10/i _1 kpc) on the source plane, we prepare new ray- 
tracing simulations with higher angular resolution of 0.1 
arcsec by narrowing the field-of-view with the fixed num- 
ber of light-rays. The field-of-view is set to be 4 x 4 
arcmin 2 , the number of light-rays is 2048 2 , and 100 real- 
izations are prepared. The resulting angular resolution 
is 4arcmin/2048 = 0.1 arcsec. There are ~ 15 light-rays 
on the circle with the radius of 10/i -1 kpc even at the 
highest redshift z s =20, which is sufficient to calculate 
the average of the magnification given by Equation (fT8|) . 

In this section, we do not include the multiple images 
formed by strong lensing in the magnification PDF, i.e., 
we use the only rays with type I. If there is a source 



including type II or III, we do not include the sample in 
our analysis. However, since almost all rays are type I, 
our results are not significantly changed. 

Figure [16] shows the magnification PDF for various 
source radius at z s = 5. The three symbols for the point 
source (black), Rs — 10/i _1 kpc (red), and Rs = 50/i -1 
kpc (blue). We find that the high magnification tail sig- 
nificantly decreases with increasing the source size, which 
demonstrate the importance of the source size effect for 
high magnification events. 

5.2. Magnification Effect on Luminosity Functions 

The lensing magnification changes the observed flux 
or luminosity by a factor of fi. Specifically, the observed 
(lensed) luminosity L b s and the unlensed luminosity L 
are related as L bs — f-L. Since the number density of the 
source is conserved with and without lensing magnifica- 
tion, we have ^ hs (L ohs )dL ohs = $(L)dL, where $^ bs ($) 
is the luminosity function (hereafter LF) of sources with 
(without) the lensing. Thus the lensed LF is written as 

$S a (L ob .)=/^^(W/*)- (19) 
J n an 

We consider LFs of the Schechter model for galaxies and 
the double power law model for quasars to demonstrate 
the impact of lensing magnifications on observed LF. 
The Schechter LF is defined as (Schechter 1976), 

= (L/L,) a exp (-L/L,) , (20) 

where a is the faint end slope, L* is the characteristic lu- 
minosity, and is the normalization. Figure [IT] shows 
the lensed Schechter LF for redshifts z s = 1 — 20. The 
solid curves are the lensed model, while the dotted curve 
is the unlensed model. Left panel is for point sources, 
while right panel is for extended sources with the radius 
of lO/i^ 1 kpc. The left (right) side in each panel is for 
the faint end slope a = — 1 (—2). We can clearly see the 
lensing effect on the Schechter LF at the bright end, the 
exponential tail is modified to the power law. This is 
because there are many faint sources, some of which are 
magnified. This feature is consistent with previous works 
(e.g., Lima et al. 2010b; Wyithe et al. 2011). As shown 
in the Figure, the lensed LF is proportional to L~ 3 for the 
bright end L 3> L* for the point sources. We can under- 
stand this asymptotic behavior analytically by inserting 
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Fig. 17. — The lensed and unlensed Schechter luminosity function for redshifts z B = 1 — 20 are denoted as solid and dashed curves, 
respectively. The left panel is for the point sources, while the right panel is for the extended sources with the radius of 10/i — x kpc (in 
comoving scale). The left (right) side in each panel is for the faint end slope a = —1 (—2). 
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Fig. 18. — The lensed quasar luminosity function. The vertical 
axis is the lensed model divided by the unlensed one. 



dP s (n)/dfi oc fi 3 (valid for fi 3> 1) to the Equation (fT9]) 
with Equation (r20|) . Then we have $ bs oc £~ b 3 s , which 
is independent of a as long as a > — 3 is satisfied. For 
the finite source size, the lensing effects become less sig- 
nificant. The lensing effect appears for bright galaxies 
L > 10L*, suggesting that wide area surveys of high red- 
shift sources are necessary to see the lensing feature in 
the LF. 

The double power law model for the quasar LF is 

$(L) = — n , (21) 

(L/LX + iL/L^ 1 

where a and (3 are the slope for bright and faint sources. 
We set a = 3.3 and f3 = 1.4 to be consistent with the 
quasar LF from the Sloan Digital Sky Survey (Croom 
et al. 2009). Figure [TBI shows the lensed LF, the ratio 



of lensed to unlensed LFs. We find that the number of 
bright sources increase by ten percents or so, which is not 
so pronounced compared with the case of the Schechter 
LF. This is simply because the magnification effect is 
more significant for steeper LFs. 

6. SUMMARY AND DISCUSSION 

We have presented high-resolution ray-tracing simula- 
tions to derive accurate PDFs of the lensing convergence, 
shear and magnification for distant objects. The resolu- 
tion of our simulations, for instance the softening length 
of the iV-body simulations of 2/i _1 kpc (comoving scale) 
and the grid size of the two-dimensional gravitational po- 
tential on the lens planes of 3/i _1 kpc (comoving scale), 
is high enough to enable direct predictions of lensing ef- 
fects on distant galaxies. In addition, we study the PDFs 
up to very high source redshifts of z s = 20. In our TV- 
body simulations, we used the different number of par- 
ticles of 256 3 ,512 3 and 1024 3 with the fixed box size of 
50/i -1 Mpc, in order to address the effect of the numerical 
resolution carefully. We have found that the sufficient 
number density of iV-body particles is necessary to re- 
duce the shot noise effect which artificially broadens the 
PDFs especially for high-redshift sources. Both numeri- 
cally and analytically we confirmed that our highest res- 
olution run with the smallest grid size of 3/i~ 1 kpc is not 
affected by the shot noise effect. We also examined the 
effects of density fluctuation beyond the simulation box 
comparing two box sizes of 50/i _1 Mpc and 100/i _1 Mpc 
and found that the effects of box size on lensing PDFs 
are not significant. 

First we have studied the convergence PDF. We have 
found that the mean convergence (k) measured in the 
source plane is not zero, as often assumed in various anal- 
ysis of cosmological lensing effects, but systematically 
has negative values. The mean convergence is found to 
be correlated well with the variance of the convergence, 
which follows (k) ~ — 2(k 2 } as expected from simple con- 
sideration. Meanwhile the mean convergence (k) mea- 
sured in the image plane is zero within small scatters 
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less than 10 -3 . Therefore, we have to use appropriate 
statistics depending on the situation. We have found 
that the modified log-normal model of Das & Ostrikcr 
(2006) reproduces our simulation results quite well, ex- 
cept at the high convergence tail where the model slightly 
underestimate the PDF. 

Next we have shown how the magnification PDFs are 
closely related with the convergence PDFs. Specifi- 
cally, we have pointed out that the simple relation fi — 
(1 — k)~ 2 , which approximates the correlation between 
convergence and magnification seen in our simulations, 
can be used to convert the convergence PDF to magnifi- 
cation PDF (or vice versa). Surprisingly, the magnifica- 
tion PDF obtained via this conversion agrees very well 
with simulation results up to very high-magnification tail 
of [a ~ 10. In light of this finding, we have presented a 
simple analytic model of the magnification PDF based on 
the convergence PDF of Das & Ostriker (2006). We have 
also presented a simple fitting formula of the shear PDF, 
which is shown to reproduce the simulation results very 
well. As explicit applications of these PDFs, in this paper 
we have discussed strong lensing probabilities and made 
quantitative predictions for the magnification effects on 
observed luminosity functions for distant sources. 

To summarize, the convergence PDF is computed by 
Equations ([8]), (j9|), and ([5]), the magnification PDF is 
computed by Equation (fTTj) via the convergence PDF, 
and the shear PDF is computed by Equations (fT2"|) and 
(fl~3|) . All these models require the variances of the con- 
vergence and shear. Under the weak lensing approxima- 
tion, we can compute them analytically as 



<« 2 > = < 7 2 > = 
r(z)r(z, z s ) 



9 

8tt' 



o H(z) 



dZ (1 + -) 2 



-co 



dkkP(k,z)[W(kr(z)0)} 2 , (22) 



where W(x) is the window function corresponding to 
the shape and size of the source or the smoothing scale 
of interest, e.g., W(x) — 2Ji(x)/x for the top-hat and 
W{x) = exp(— x 2 /2) for the Gaussian window function. 

Our numerical simulations consider dark matter par- 
ticles only, while the baryonic effects become important 
in the galactic scale. For instance, Hilbert et al. (2007) 
investigated the effects of stellar masses in dark halos on 
the ray-tracing simulation through the Millennium sim- 
ulation. They used semi-analytical galaxy models from 
the halo merging history to investigate the baryonic ef- 
fects, and showed that the baryon increases the magni- 
fication PDF by a few ten percents at /x = 100 for the 
source redshift of z s — 2.1, and enhance the strong lens- 
ing probability by an order of magnitude at z s > 1. We 
are planning to include the baryonic components by sim- 
ply placing galaxies (bulge and disk components) in dark 
halos. Such simulations should also allow us to study cor- 
relations between the magnification map and the distri- 
bution of foreground objects such as galaxies and clus- 
ters (Takahashi et al. in preparation). Similarly, the 



dark matter substructures would affect the small-scale 
power spectrum and the weak lensing observables. In 
fact, Hangan et al. (2005) showed that the subhaloes en- 
hance the power spectrum at k > 100/i(Mpc) _1 . Hence, 
the unresolved subhaloes in our simulation could affect 
our prediction. 

It is known that the statistical property of the random- 
Gaussian field is fully characterized by the power spec- 
trum. However, the characterization of the non-Gaussian 
fields, such as the lensing field studied in this paper, re- 
quire information on higher-order correlation in addition 
to the power spectrum. In fact, it has been shown that 
the power spectrum contains little information in non- 
Gaussian regime (e.g., Rimes & Hamilton 2005). There- 
fore, the PDFs studied in this paper are expected to 
contain useful additional information on the statistical 
property of the lensing map. Recently, Neyrinck et al. 
(2009) showed that gaussianizing the one-point distribu- 
tion function of the matter density fluctuation by using 
a modified log-normal transformation, S — > ln(l + 5), in- 
crease the signal-to- noise ratio (S/N) of the transformed 
power spectrum. The result implies that such gaussian- 
ization provides a means of extracting non-Gaussian in- 
formation and recovering the information content of the 
power spectrum. Similarly, Seo et al. (2011) showed that 
the cosmological information in the convergence power 
spectrum is recovered by using a modified logarithmic 
transform of the convergence field. The accurate PDFs 
presented in the paper might be useful in this regard. 

Our simulation results, including two-dimensional 
maps of the convergence, shear, magnification, and 
PDFs of these quantities, are publicly availa ble at 
http://cosmo.phys.hirosaki-u.ac.jp/takahasi/raytracing/ 
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APPENDIX 
BOX SIZE EFFECT 



Throughout this paper, we have shown the results for the simulation box 50h 1 Mpc on a side (Model S256, S512 
and S1024, see Table [1]). However, the density fluctuation larger than the box size may affect our results. In order 





Fig. 19. — Comparison of two simulations with a different box size, the small box L = 50h Mpc (model S256) and the large box 
L = 100/i _1 Mpc (model L512). These two simulations have the same mass and force resolutions (see Table [TJ. Top left panel shows 
the root mean square of the convergence, the dots are the results for S256 and the circles are for L512. Top right panel plots the strong 
lensing probability. Bottom panels display the magnification PDF (bottom left), the convergence PDF (bottom middle), and the shear 
PDF (bottom right). The black symbols are for S256, while the red symbols are for L512. These result indicate that the box size effect is 
negligible. 



to check such potential box size effect, we compare our simulation results with those of the larger simulation box 
of lOO/i -1 Mpc (Model L512). We compare the results for the model S256 (smaller box) and L512 (larger box) in 
Figure [TH1 Note that the two models have the same mass and force resolutions. We check the root mean square of the 
convergence (top left), strong lensing probability (top right), and PDFs of magnification (bottom left), convergence 
(bottom middle) and shear (bottom right). We find that two results agree well with each other. A possible exception 
is the strong lensing probability and PDFs at the high magnification and convergence tail, which appear to be slightly 
enhanced in the larger box size simulation. The differences is presumably because more massive halos are formed for 
larger box simulations. In either case, our results here confirm that the box size effect is insignificant. 
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